Symmetry, Integrability and Geometry: Methods and Applications 



SIGMA 8 (2012), 033, 13 pages 



A Top-Down Account of Linear Canonical Transforms* 

Kurt Bernardo WOLF 

Instituto de Ciencias Fisicas, Universidad Nacional Autonoma de Mexico, 

Av. Universidad s/n, Cuernavaca, Mor. 62210, Mexico 

E-mail: bwolf@fis.unam.mx 

URL: http : //www. f is .unam.mx/~bwolf / 

Received April 24, 2012, in final form June 01, 2012; Published online June 06, 2012 
http://dx.doi.org/10.3842/SIGMA.2012.033 

Abstract. We contend that what are called Linear Canonical Transforms (LCTs) should 
be seen as a part of the theory of unitary irreducible representations of the '2+1' Lorentz 
group. The integral kernel representation found by Collins, Moshinsky and Quesne, and the 
radial and hyperbolic LCTs introduced thereafter, belong to the discrete and continuous 
representation series of the Lorentz group in its parabolic subgroup reduction. The reduction 
by the elliptic and hyperbolic subgroups can also be considered to yield LCTs that act on 
functions, discrete or continuous in other Hilbert spaces. We gather the summation and 
integration kernels reported by Basu and Wolf when studiying all discrete, continuous, and 
mixed representations of the linear group of 2 x 2 real matrices. We add some comments on 
why all should be considered canonical. 
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1 Introduction 

Linear canonical transforms (LCTs) have been developed as the diffraction integral kernel for 
generic paraxial optical systems by Stuart A. Collins [7], and also defined as the group of 
unitary integral transforms that preserves the basic Heisenberg uncertainty relation of quan- 
tum mechanics in D = 1 or higher dimensions by Marcos Moshinsky and Christiane Quesne 
[16, 17, 23]. Further, LCTs can be seen as the group actions generated by the Lie algebra of 
quadratic Hamiltonian operators [25]. Also, radial [18, 26] and hyperbolic [27] canonical trans- 
forms have been defined after separation of variables in D > 2 dimensions. More recently there 
has been interest in canonical transformations that are represented by matrices, infinite, semi- 
infinite or finite-dimensional - the latter due to the finite capacity of measuring and storing 
devices, but with serious concomitant difficulties. 

In both the optical Lagrangian formulation [7] or in the quantum-mechanical Hamiltonian 
approach [16], one arrives at a parametrization of all possible linear transformations through 
a 2D x 2D real symplectic matrix M E Sp(2D, R). For reference we gather here the basic 
D = 1-dimensional formulas. We consider the matrix 



as the presentation of the group Sp(2, R), which is equal to the groups of all real 2x2 matrices 
of unit determinant SL(2, R), and is isomorphic through a complex similarity transformation to 
the group SU(1, 1) of 2 x 2 pseudo-unitary matrices. This group of 2 x 2 matrices in fact covers 




ad — be = 1, 



(1) 
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twice, i.e., is 2:1 homomorphic to the '2+l'-Lorentz group of 3 x 3 pseudo-orthogonal matrices 
L(±M) € S0(2,l), 

, ; (| , 2 {a 2 -b 2 -c 2 +d 2 ) bd-ac \(a 2 -b 2 +c 2 -d 2 f 

^— > I cd—ab ad+bc —cd—ab I . (2) 



c d 



J 1 {a 2 +b 2 -c 2 -d 2 ) -bd-ac \{a 2 +b 2 +c 2 +d 2 ) 



These accidental iso- and homomorphisms between symplectic and relativity groups have been 
exploited to describe and relate various physical and optical models [11, 15, 19, 31]. 

However, in spite of its apparent simplicity, the group of 2 x 2 matrices (1) has an infinite cover 
group Sp(2, R). It is the double cover - the metaplectic group Mp(2, R) - which is represented 

by the well-known integral form of linear canonical transforms C M = C[ a , ) of the usual Hilbert 



space of functions f(x) £ L (R), 

/mW = (C M f)(x) = / dx'C M (x, x')f(x'), 



c d 

2 



JR 

C M (x, x') := 1— exp ( ^- [dx 2 - 2xx + ax 2 ) ) , (3) 
y2mb \2o 

where the phase of the prefactor is taken to be 

1 1 



\/2^i6 y/2ir\b\ 



exp (— ij7rsign&) . 



Unitarity is evident in that C M — i (x, x') — Cm (^'? x)* . 

Phases are quite delicate here; Collins did not consider the group property of the trans- 
forms [7]; yet Moshinsky and Quesne [16] realized that C Ml C M2 = ±C Mi m 2 > a sign depending 
on the signs of the bi, 62 and 612 elements in a rather complicated way [29, Chapter 9]. The 
double-cover issue can be simplified by observing that for F = , the kernel (3) is that 

of the Fourier transform J 7 , but for a phase: C F = e~ m I^T . Since J 74 = 1 , only C\ will return 
the cycle to the unit 1 . In the limit 6 — ^ from the lower complex half-plane, the kernel (2) 
becomes a Dirac delta, 

v ^ / i\ exp(icx 2 /2a) / , x\ 

hm C M (x,x ) = n — -Six . (4) 



&->o ^/a V a 

We may skip further detailed consideration of the D = 1 integral linear canonical transforms 
and their properties, which are mostly standard knowledge, and from whose D = 2 case one can 
build the radial and hyperbolic LCTs in a bottom-up construction. Rather, the purpose of this 
review is to give a top-down panorama of LCTs. 

In Section 2 we return to the Lie algebra so(2, 1) of the Lorentz group S0(2, 1) realized 
by second-order differential operators corresponding to harmonic and repusive oscillators, with 
a singular centrifugal or centripetal potential, and the generator of scaling. Linear combinations 
yield free propagation and a square-radius coordinate. In Section 3 we list the eigenfunctions 
and spectra of those operators for the centrifugal case that fall into the Bargmann T>k discrete 
series of representations. There are three subgroup orbits that are examined in three subsections: 
elliptic, parabolic, and hyperbolic. Section 4 follows the same structure for the centripetal case 
that fall into the Bargmann Cf continuous series. This is what we deem to be six faces that 
Sp(2, R) linear canonical transforms can show in various Hilbert spaces. The resulting matrix 
and integral kernels were computed in 1981 by Basu and Wolf [6], so the results are not new, but 
are gathered here for the first time as proper, unitary LCTs. The concluding Section 5 offers 
some reasons to call all representations canonical, and some comments on finite-dimensional 
approximations to these transforms. 
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2 The Lorentz algebra 

The generic form of the second-order differential operator realization of the Lorentz Lie algebra 
so(2, 1) can be obtained from the oscillator algebra [17], adding a term 7/r 2 to the second 
derivative, whose interpretation is that of a centrifugal potential as 7 = m 2 — \ ~> —\ for 
solutions of angular momentum m G Z. Separation into hyperbolic coordinates yields 7 < — \ 
when added by a centripetal force. The formal operators are associated with Hamiltonians, and 
are generators of Sp(2, R) subgroups as follows: 

n == \ (-£ +}< +<>) ■ «p0«n - (£# T4:) • 

where the 'harmonic oscillator' Jq generates the compact {elliptic) S0(2) subgroup (and its co- 
vers), while the 'repulsive oscillator' and 'scaling operator' generate equivalent noncompact 
(hyperbolic) subgroups S0(1, 1). Their commutation relations are 

[J7, j]] = -ui [ji jgi = iJ7, [ji J7] = ui (6) 

Also relevant are the linear combinations 

n ■= ^ +JJ = \ (~ 2 + j) , exp(iWj) » (J "*) , 

J2:=^-J7=l r 2 , exp(icJ2) m. ^ J), (7) 

where the 'free system' and the 'square-radius' JZ generate the Euclidean (parabolic) sub- 
groups ISO(l), which are equivalent under the Fourier transform. Note carefully that these are 
not 'raising and lowering' operators for discrete eigenvector bases. Instead and particularly, 
J2 determines the diagonal radial position operator with respect to which we refer the eigen- 
functions of all other operators, as done below. These eigenfunctions depend crucially on the 
value of 7 in (5), as given by the Casimir eigenvalue 

C:= (Jl? + (J2f-(J2? = KU 

K = -b + re=- k ^-k), fc = |(i±^T + |Y (8) 

where k is the Bargmann index [2] that determines (up to parities) the essentially self-adjoint 
irreducible representations of the algebra in a Hilbert space £ 2 (R + ) with measure dr - which 
are unexpectedly imbricate. 

3 The T>k canonical transforms 

We consider first the case when the coefficient 7 is of centrifugal origin in D = 2 using polar 
coordinates x = rcos9 and y = rsin#, r G R + , 9 modulo 2n. We want the measure to be dr, so 
after similarity transformation, v / rV 2 /y / r = d 2 + r~ 2 (\ + <9 2 ), with <9 2 1— > —m 2 for the subspaces 
of angular momentum m G Z so that k = ^(\rn\ + 1) G {|, 1, |, . . .}. This is the Bargmann 
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discrete representation series T>^ [2]. 1 The eigenfunctions and spectra of the five operators (5) 
and (7), are [6] 



( r ) = J — r 2fc-l/2 -r 2 /2 L (2fc-l) f 

rn( ' y (2fc-n-l)! ' 1 j 

V2(2fc+ra-l)! 2fc _ 1/2 2/2 /-n_ 2 



r(2fc)Vn! V 2fc 
■$JM = e™ k ^J 2k -i{pr), p e R+, 



! r ^-i/z e -r-/^ lFl ( "-.^ ) 5 ,„ ... /, + // . „ a Z. 



-**(r) = 5(p-r), p€R+, 



lAfc^ _ e^ (2fc+M)/2 2^r(/t+i/i) 1 



W) = w^S /= ^ ^M iM _ 1/2 (-iH), 



= r (2fc )^ r e lFl { 2k >- ir ) ' ^ e 

2 <^(r) = -Lr- 1 /^, /, € R, (9) 

where are the Laguerre polynomials, J u is the Bessel function of the first kind, and M\ „ 
is one of the Whittaker functions. In Dirac notation, x $k(r) = _(£;, r| A)x, with X G 
{0, 1, 2, +, — } indicating the eigenkets of Jj, with eigenvalues A G £x in the spectra of (9). 

Next, we should obtain the transformation of the functions in (9) under the generic LCT 
C M = To this end for each generator in the lists (5) and (7) we decompose C M into 

a right-ordered product of the subgroup with a diagonal-plus-phase transformation as (4), 

C( a c b d )=c(t ^Je^M)- (10) 

where a', d and a will be algebraic and trigonometric functions of a, b, c, d, respecting ad— be = 1. 
The right factor applied to $|(r) will multiply it by the phase e lAo , while the left factor will 
be given by (4), with the sole precaution of taking = (sign a) 2k /^/\a\ stemming from the 

radial reduction of the D = 2 case with angular momentum. 



3.1 Elliptic basis: lower-bound discrete X^-LCTs 

We consider first the eigenbasis of the compact generator Jq, $^(r), in the Bargmann 
representation series, whose eigenvalues are lower-bound and discrete, X = m = k + n, n *E Zq . 
The decomposition (10) yields a' = \/ a 2 + b 2 , a'c' = ac + bd and e ia = (a—ib)/(a+ib). Then, 

„2 ac+bd 
2(a 2 +b- 

a+ibj (a 2 + 6 2 )V4 * m \tffi+p 
From this we can find the representation matrices of the group 



c( : 5) v.w - ( ^ m 6 T. ^ V. I I ■ (ID 



^ A ,(M)= / x^Al^r^dr^rlCMl^AV, 
J o 



The same were named by I.M. Gel'fand and M.A. Naimark to be the complementary representation se- 
ries [8, 21]. Note that k o 1 — k correspond to the same k; the interval < k < 1 is exceptional in that 
the operators have a one-parameter family of self-adjoint extensions [30], and also harbor the exceptional, or 
supplementary irreducible representation series [2, 8, 21]. We shall be working within the Friedrichs extension. 



A Top-Down Account of Linear Canonical Transforms 



5 



by straightforward integration of the confluent hypergeometric functions that appear in the 
tables of Gradshtem and Ryzhik [10]. In the SO(2) basis, this is a matrix with rows and 
columns numbered by m = k + n, n G Zq, that was found by Bargmann [2], 

°<m'(: 5) : = (°<> C (c =J o °°dr% k m (ryc( a c b d )% k m ,(r) (12) 

2 2k T(m+m') 
^T(k+m)T(l-k+m) F (k+m')T ( 1-k+m') 
x [(d-a)-i(6+c)] m - fc [(a-d)-i(6+c)] m '" fc [(a+rf)+i(6-c)]- m - m ' 

/k-m, k-m' a 2 +b 2 +c 2 +d 2 +2 \ 
X 21 V l-m-m> ; a 2 +6 2 +c 2 +d 2 -2 y ' [ ' 

Then, when f = {/ n }^=o ^ s a vec t° r hi the Hilbert space of square-summable sequences £ 2 (Zq), 
the Sp(2, R) action 

oo 

f M EE C M : f = °D fc f, f M , n EE (f M ) n = £ H + n,k+ni(M)f n *, (14) 

n'=0 

is an LCT which is unitary in the Hilbert space of sequences £ 2 (Zq). 

In fact, (13) provides a unitary summation transform kernel for every value k > 0, not 
necessarily stemming from integer angular momentum fn, where k = ^(\m\ + 1) G {|, 1, |, . . .} 
are single- and double-covers of the SO(2, 1) group, i.e., single covers of Sp(2, R). When k is 
quarter-integer (in particular k = \ and |, to be revisited below), we have representations of 
the metaplectic group Mp(2, R), and various higher covers for fractional k's. 

In the Bargmann representations, the generator Jq has an equally-spaced spectrum 
{ m }m=k wmcn is bound from below by k > 0. There is also a paired series of representations T>^ 
where that spectrum is upper-bound by —k, i.e., m = —k — n, n £ Zq". This stems from the outer 
algebra automorphism Jq — Jq (reversing the sign of the spectrum), Jj — J^, «-> J^ j 
and so J^ — J±- This generates an outer group automorphism whose representations in any 
subgroup reduction yields the T>^ matrices or integral kernels 

Xr,k,(-)fa b\ ^X n kX+) ( a -b 



where J^ has spectrum {A} G and parity a € {+, — } under this automorphism. In fact, 
the matrices in (15) are related by a similarity transformation with ^q ^j, a non-symplectic 
matrix which describes reflection in geometric optics [28, Chapter 4]. 

We have spoken of T>^ for k > 0. What about k < 0? The Bargmann index k acts as a 
lower bound for the values of m in T>^ , and as an upper bound in T>~^ in the one-step recursion 
relations obtained with raising and lowering operators on the eigenvectors of Jj; but they are 
one-way barriers. At k = the m = ±k lines cross, so that when k < is integer or half- 
integer, the m's have a lower bound at negative k and an upper bound at positive —k. Raising 
and lowering the m's in this range will yield a (2/c+l)-dimensional matrix, which is a faithful, 
although non-unitary irreducible representation of Sp(2, R). Of course, a well-known theorem 
(see [9]) states that non-compact groups do not possess finite-dimensional unitary irreducible 
representations. 

3.2 Parabolic basis: continuous radial X^-LCTs 

The eigenfunctions + <& k {r) and ~<& k {r) in (9) of the parabolic generators Jj_ and JZ in (7) form 
generalized bases for £ 2 (R + ). Consider thus the Sp(2, R) representation given as an integral 
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kernel by 



DlA: b d ) ■= MM: 5 r*£ = MM: s c ! ~o + *M m 



P,p'\c d) • \ P'M C d I ^ p' I \ ^P'Mc ci/Ml / P 



(-** c 



* 1/6/(0 "f + ^ 



e-hr* (. df?+ap! 2 \ r (pp'\ 

■ 'pf/exp 1 — J 2fc _x — (17) 



b v rr 2 & / V & 

_ 2{ Pp i? k - 1 ' 2 pxn ( df?-2 P f/+a^ \ (2k — \ 2ipp'\ 

This we recognize as the radial canonical transform kernel [18, 26] for angular momentum k = 
\(jn + 1), which can be extended to k > 0, that acts unitarily on the functions f = {/(r)} rg R+ 
in the Hilbert space £ 2 (R + ) with measure dr. As we had in (14), now 

f M = C M : f = -D fc (M)f, f M (r) = (f M )(r) = / dr , D* r ,(M)/(r / ). (19) 

JO 

We regain the classical D = 1 LCT kernel in (3) as the direct sum of the representations 
k = 7 and |, for functions /(x) whose domain is extended to x £ R through writing them 

as the sum f(x) = fe l ^\x) + fo^ 4 \x) with even and odd parity summands f e (—r) := f e (r) 
and fo(—r) := —foir) respectively. For 2k— 1 = =pl, the integral kernels contain J_]/2(z) = 
^2/ttz cos z and J + i/ 2 (.z) = \/2/irz sin z, z = pp'/b, with phases i 2fc = e 17r//4 and ie 17r//4 . Their 
sum thus yields e 17r / 4 e _1:r:!; / b i n the oscillating Gaussian of the original LCT kernel (3), 

C M (x,x') =Diy 4 \r,r') + D^\r,r'). (20) 

Had we chosen the eigenbasis of the free Hamiltonian Jj_ instead of the square-position Jl , 

the Bessel function + & k (r) in (9) we would have the Hankel transform of (18) by F = JJ, 
so 

+d p,p>{: 5) : = ( + <^ c (: 5) + ^) = "^V(-6 r)- (21) 

3.3 Hyperbolic basis: the face of X^-LCTs 

There remain the eigenbases of the two equivalent noncompact operators in the list (5): the 
repulsive oscillator Hamiltonian and the scaling generator ■ The latter is the simpler of 
the two because its eigenfunctions, 2 & k l (r) in (9), are the Mellin transform kernel with p G R 
and k > 0. Since we know the LCT action C M in the parabolic basis, (16)-(19), we apply it to 
these eigenfunctions, 



(C M : 2 ^)(r) = ^[^'^%{: b d y- 1/2+itl 

e~ ink T{k+ip)r 2k - l l 2 e ldr2 / 2h ^ ( k+ip -ir 2 



2k-ip^ Y(2k) b 2k (-\a/b) k+ '^ \ 2k ' 2ab 



& o r;^r • (22) 



Needless to say, the joint phases of —ia/b are needed so as not fall into multivaluation problems. 
This result, reported in [6], was calculated following the general method to find Mellin transforms 
of hypergeometric functions due to Majumdar and Basu [4]. 

Now it is only necessary to perform the Mellin transform of (22), to obtain 



2nfc (a b \ . (Zfak nl a b \ 2,f,fc 

U p,p'\c d — [ > L c d V 
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r(fe-i/x)r(fc+i//) 

2irT(2k) 



fey y 2k ad 

One should note that the complex power functions are evaluated along the imaginary axis, in 
the principal sheet where the cut is chosen along the negative real half-axis. 

Let us call this the P^-hyperbolic basis (not to confuse it later on with hyperbolic LCTs), 
because it lies the T>t representation of Sp(2, R). The corresponding LCT with the integral 
kernel (23) that transform functions f = {/(//)}u 6 r E £ 2 (R) unitarily, are 

/oo 
d// 2 Z^,(M)JV). 
-oo 

Since and Jg are related by similarity through Q ~^ J (the square root of the Fourier 
transform), we obtain an equivalent integral LCT, 

to* , (° ^ = to* ± f a+ J +C+ ^ - a+b ~ C+d ] . (24) 
^ dj ^ 2 \-a-b+c+d a-b-c+d J y ' 

In this section we have thus shown three faces of LCTs: the summation kernel (13) for vectors 
{/n} ngZ +, the 'radial' canonical transform integral kernel (18) for functions {/(r)}reR+ and its 
Hankel transform (21), and the (apparently unknown) P^-hyperbolic transforms (23) and (24) 
for functions {/(/x)}ugR. These various functions are in fact the coordinates of the same abstract 
vector f as f n = o(k, k+n\f), f(r) = r|f), f(p) = i(k, /x|f), etc. In this sense, all faces of 
the Hilbert space vector f are related, and so are their LCTs. 



4 The C e s canonical transforms 

Let us now consider the range 7 < —\ of the so(2,l) generators in (5), corresponding to cen- 
tripetal potentials, which also stem from the separation of D = 2 coordinates in two disjoint 
patches of hyperbolic coordinates. One patch (indicated by a = +1) is x = pcosh£, y = psinh£, 
and the other (a = —1) is x = /)sinh£, y = /9cosh£, for p,( 6 R. Fourier expansion in Q pro- 
vides d 2 1— > —(, , ( 6 R that yields a centripetal term in the generators (5). In contrast to 
radial coordinates, which follow the subgroup reduction Sp(4, R) D SO(2) ®Sp(2, R), hyperbolic 
coordinates conform to the reduction Sp(4, R) D 0(1, 1) (g) Sp(2, R), where 0(1, 1) contains the 

discrete reflections ^ J J that interchanges the two values of a, and that allows the 

reduction by even or odd functions in p, so the range of this hyperbolic radius is reduced to 
r = \p\ £ [0, 00) and a = sign p. The Hilbert spaces of functions for the range of representations 
7 < consists of functions _(r, a\f) = f(r,a) = f a (r) whose inner product we can represent 
with two-vector notation 

i{r) = ( f f +1 ) (r) e ^(R + ), (f,g) := £ I™ drU(r)*g a (r), 

taking care to note that a is an index of the basis functions |r, er)_ stemming from the two 
hyperbolic coordinate patches. When 7 = — i, the Bargmann index in (8) is k = \\ for 7 < —4, 
k becomes complex, 

k = 5 + is, seR, k = \ + s 2 >\ 
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and determines these representations to belong to the Bargmann continuous nonexceptional se- 
ries Cg, where e G {0, ??} distinguishes between vector and spinor (two-fold cover) representations 
of SO(2, 1), and in the SO(2) reduction the index m = e modi is unrestricted 2 . 
The (reduced) list of eigenfunctions of the generators J\ in (5) is now (cf. (9)), 

o*e,*( T ) = MM ((-l) m ~ e V^(k-m)T(l-k-m)W m ^ 1/2 (r 2 ) 
m y ' vr^ V 2T(k+m)T(l-k+m)W„ mik _ 1/2 (r 2 ) 

m+e G Z, k = \ + is, go(k) = coshes, ffi/2(&) = sinh7rs 



*?{r)=r\ M forp>0, and " for p < 0; (25) 



**;*(r) = --L r- 1 ^, ^ G R, t G {+1, -1}; 



where W\ „ is a Whittaker function. The two components of the eigenfunction are distinguished 
by the non-symplectic matrix G 0(1,1), which on the algebra generators (5) and (7) 

can be interpreted as 2 x 2 matrix operators Jq = ( _ j 7 J , = f J q _°p ^ , J2 
J ± - -Jg 



4.1 Elliptic basis: discrete C^-LCTs 



The action of the LCT operator C M on the discrete eigenfunction basis of Jq, (r), for 

m — e G Z will be that of an infinite matrix C^J = ||C^, / 1| . The transformation of this basis 
under Sp(2, R) can be again factorized into the right-ordered product in (10), so that as in (11), 
we now have 



(r) - I 6XP v 1<Tr2 ^I/ 0^, fc / ; 1 



where the component sign <r only appears in the Gaussian exponent. This relation shows that, 
in the same way as harmonic oscillator functions reproduce under the LCTs (3), and Bessel 
functions under radial LCTs (16), the Whittaker functions in (26) do likewise under the Cf 
continuous-series LCTs. 

But the following step of finding the LCT matrix elements as in (12), results in a sum of 
integrals of two Whittaker functions and an oscillating Gaussian, which the authors [6] could 
not solve. As in (13) this result was obtained before with the traditional m-shift operators by 
Bargmann [2], 




formV 2 2m l T(k+m)T(l-k+m) [{a-d)+i(b+c)] m - m 
T(k+m)T(l-k+m) [(a+d)+i(b-c)] m+m ' 

h 771 , 1 k 771 1/9 ,9 9 ,9 ^- 

;-\ a 2 +6 2 +c 2 +d 2 -2 

1+m—m' 



for m<m' 



("I) 



m —m 



2 2m r(k+m')T(l-k+m') [(a-d)-i(b+c)] m '- 
"mTV T(k+m)T(l-k+m) [(a+d)+i(6-c)] m ' H 



2 I.M. Gel'fand and M.A. Naimark called this the principal representation series [8, 21]. 
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X 2*1 ™? \ k m ;-l(a2 +62+c 2 +d 2_ 2) 

\ 1+m — m 

As in (14), but now for infinite vectors of discrete components, f = {f m }m-eeZi the LCT is 

oo 

m'=— oo 

and is unitary in the Hilbert space of sequences i 2 (Z). 
4.2 Parabolic basis: radial CJ-LCTs 

The action of the canonical transform operators C M on functions _(<t, r|/} = /<r(?") £ wm 
be represented in the continuous series Cf by a 2 x 2 matrix of integral kernels 



f M (r) 



r, r') 



f+i 
f-i 



dr' C £ ^ k (r,r')f(r'), 



r ') Cm;+1 ,-l( r ' r ') 
" C ?-l,+l( r i r ') ~ ^-1,-1 ( r ' r ') , 



(27) 



Here we resort to the reduction Sp(4, R) D 0(1,1) (8) Sp(2, R) as presented in [27], with the 
precisions made in [6, equations (2.15)] regarding ranges and phases. For the non-exceptional 
continuous series of representations C e , where k = k(l— k) >\ and 



£ = : he 

c - I . J, 



2 ' 



, fc— 2 = is, s > 0, 
-1, fc— \ = is, s > 0, 



(28) 



where g e (k) is given in (25). 



The integral kernel elements of C^* are then expressed as 



( C^J ) a y(r,r') = G M;CTi(T /(r,r / )^' (T ,(-rr76), 
_ , A -v/rr 7 f .dor 2 + ao'r' 2 \ 

GM -' (r ' r):= 2 7 TH exp V ~J' 

= M*!£_i(C) = heH^ +l (-C) = H$£{C), 
H+1,-i(0 ■= 4(-signC) 5£ (A;)K 2is (|C|) 

= ^fl!} t+1 (C) = kfl?f,_i(-C) = M^t-iCO, 



(29) 
(30) 

(31) 
(32) 



where and are the Hankel functions of the first and second kind valued above and 
below the branch cut, K v is the Macdonald function, while g e (k) and h £ {k) are given in (28). 

The 2 x 2 kernel LCT in (27), Cf M ' aa ,(r, r'), r,r' G R+, can be written also as (p,p'), 
with p := or G R. In that case, the kernel C e s has been expressed [6] in terms of confluent 
hypergeometric functions, 



Cm (p,fJ) 



(signb) 2 *h £ 1+sisap ' )/2 g £ (k) 



n\b\ 

r(i-2jfc) 



exp i 



dop — 2-qpp' + ao'p' 



I J 2 



2b 



PP 



PP_ 
2b 



2k-l 



2k-\ 2\pp' 
r\b 



2 

Ak-l 



+ [k^ 1-k] 



(33) 
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where a = sign p, a' = sign p'\ rj = 1 for a = a' and rj = — i for a ^ a'. In particular, for 6 = 0, 



ek (sign d) ( Cp^\ 

Hiti-Cm' (p,p) = 7 = r -exp isignp— S(p' - p/\a\). 

6-+o ^/|a| V 2a / 

We are not addressing the Bargmann exceptional representation series, for which in (28) one 
has k = 9 + s, with < s < \ overlapping the exceptional interval, and for which (29)-(32) 
appears to be a valid integral transform. Treatment of the exceptional continuous series can be 
seen in [5] . 



4.3 Hyperbolic basis and their C^-LCTs 

The eigenfunctions 2 &r'^ l (r) of in the continuous series representations form a generalized 
basis for a Hilbert space C?,(R), and are given in (25) with r G {+1,-1}, p G R and r G R + . 
They can provide an LCT form with a 2 x 2 matrix of integral kernels that are the double 
Mellin transforms of the LCT kernel in the position basis given by (29)-(32) or (33). Using the 
technique of Basu and Majumdar on the latter, [6] reports 



2s-iE,k fab 



-signfr) i£ g £ (fc) 
2vr 



TT'h 



rh f 



where 



ctk + — - + r'/3 k + , 

Cik Pk 

+ [ h £ a\-k + + r'Px-k + J Ti_ k 
Ctl-k PX-kJ 



r(l - 2k) Y(k - i/j) T(k + ip') fk-ip, k+ip' _ 1 
\a\ k ~^'\b\^-»')\d\ k -'w 21 V 2/c 'ad 



Tk 



ak := exp (i^7r[(/c + \p) sign ab + (k — ip) sign bd]) , 
f3k := exp (i^7r[— (/c + i//) sign a6 + (k — ip) sign bd]) , 

and where <? £ (fc) and h £ are given in (28). 

Note that the two components of this C\{R) are distinct from the two components of the 
Hilbert space C^iR^) in the parabolic case of the previous subsection. In this hyperbolic C £ s 
representation, LCTs are maps of two-component functions by the integral 



f = C ■ f — 2 P £ . fc f 

i M — L- M .i — i , 

/m;t(A*) = (f M )r(A t ) = 



,'£{-1,1} X 

and is unitary in that Hilbert space. 



/oo 



5 Concluding remarks 

The recompilation of the six forms that LCTs have in the elliptic, parabolic and hyperbolic 
subgroup bases of Sp(2, R) in the T>k and C e s representation series, has been made for the purpose 
of placing the better-known LCT forms, i.e., the 'linear' (3), the 'radial' (17), and the (lesser 
known) 'hyperbolic' (29)-(32), in the general context of group representation theory. We still 
have to justify the appellative of canonical though, because it is a term associated with the 
preservation of area elements in classical phase space, the conservation of energy in paraxial 
wave optics, and of uncertainty in quantum mechanics. 
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Classically, the 2x2 symplectic matrix M = (® in (1) acts on the phase-space coordinates, 

written (r, p) T as if they were a 2-vector. The correspondence (2) carries its action as a (2+1)- 
Lorentz transformation of the 3- vector of phase space functions 

(6),6,£2) T := (l(p 2 +l/r 2 +r 2 ), \{p 2 +7 / V 2 -r 2 ) , \rp) T , (34) 

which close under Poisson brackets into the same algebra, with i[-, •] i— > {•,•}, as the operators 
{J7}f =0 in (6), 

{&,&} = -&, {6,6} = £1, {Co,6} = 6- 

We may consider these relations as the Berezin brackets defining the Lorentz algebra. The 
square length of the vector (34) is £ 2 + — £q = — j7 [instead of — ^7 + ^ in the operator 
case (8)], and this surface is conserved under Lorentz transformations. For centrifugal 7 > or 
centripetal 7 < these surfaces are respectively a two-sheeted or one-sheeted hyperboloid, and 
a cone for 7 = 0; these surfaces are also symplectic manifolds. Since classically r = iy^^o - £1) 
and p = a linear transformation in the £'s is nonlinear in (r,p) T , Moshinsky termed these 

nonlinear canonical transformations [14, 18]. 

We count the Sp(2, R) canonical transforms collected here as different faces of LCTs acting on 
different function Hilbert spaces: £ 2 (Zq) and £ 2 (Z) in the elliptic basis of T>~£ and C e s representa- 
tions respectively; similarly, £ 2 (R + ) and .C 2 ,!^) m the parabolic basis, and C 2 (R) and /^(R) in 
the hyperbolic basis. Of course, the best known ones are those in the eigenbasis of Jl = \r 2 un- 
derstood as the position observable, which include the 'radial' LCTs in the particularly the 
'linear' transform kernels written by Collins [7] and Moshinsky et al. [17], which are genuinely fit 
to describe the signals or wavefunctions that traverse a paraxial optical or quantum mechanical 
system. For the same parabolic subgroup, the 'hyperbolic' LCTs in the C e s representations have 
found no proper application, save their use in [6] to find all SL(2, R) group representations in all 
subgroup and mixed bases reported there. Widely used in other contexts, the Bargmann results 
of 1947 [2] that provide the matrix representations of the Lorentz group, uses the elliptic sub- 
group to provide the row indices. This has been used in [1] to propose a covariant discretization 
for axis-symmetric D = 2 systems, where the limit of the discrete to the continuous radial model 
is pointedly addressed. Among the six LCT forms presented above, similar limits from elliptic 
or hyperbolic to parabolic subgroup bases should occur in all representations. Representations 
referred to the hyperbolic subgroup were investigated by Mukunda and Radhakrishnan [20] but 
have also failed to be associated with some model of physical system. 

Uncertainty relations are preserved under all LCTs, but they stem from so(2, 1) - not from the 
Heisenberg-Weyl algebra - and involve the coordinate \r 2 > 0. Recall that the mean of a self- 
adjoint operator J in a wavefield if) is := (ifi, Jtp), with the inner product appropriate to its 
Hilbert space; its dispersion is then A^(J) := ||(J— J^)ip\\ 2 - Note also that from the commutator 

between J_ = \r 2 and its (Fourier) 7^ transform operator J + , is [J_, J + ] = 2U2, and it 
follows that 

A^(J_)A^(J + ) = A^(J_)A^(J_) > J 2 ^)| 2 , 

where ip = "qMV 1 - In the special case of the original LCT (3) identified by (20), this can 

be written in the form of a Robertson uncertainty relation {r 2 )^{p 2 )^, > j((?T+P r ))^ [24]- 

We have not addressed in any depth the exceptional interval < k < 1, i.e., — \ < 7 < | 
characterizing weak centripetal and centrifugal potentials, which Bargmann treated with a non- 
local measure [2, 5]. Also there, the discrete and exceptional series overlap, and its generators J\ 
have one-parameter families of self-adjoint extensions [30], with the result that their spectra are 
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generally not equally spaced, save for the Friedricks extension in the T>k representation series. 
The various limits to the point k = \ are particularly troublesome [6]. 

Perhaps most important is the computer use of LCTs to model physical systems, pursued 
since the 1990s, in a quest that is running into hundreds of references [3, 12, 13, 22] where 
finite data sets or pixellated images are subject to LCT matrices of N x N. Of course, we know 
that noncompact groups such as Sp(2, R) cannot have finite-dimensional unitary representations. 
So, approximations have to be made, either through sampling functions and kernels, or using 
the finite-dimensional k < representations of the T>k series mentioned in Section 2; in both 
strategies the matrices will not be unitary, but only in the second will they faithfully represent 
the LCT group Sp(2, R). Perhaps the best approximation strategy is to use a truncated elliptic 
basis { *&(r)}£JJ -1 of functions that are radial oscillator modes; the matrices would by neither 
unitary nor would faithfully represent the group; but if the signal contains mostly low energy 
modes it could be fairly approximated by the first N modes, and with an inherent ./V — > oo limit 
to a representation that is both unitary and faithful. 
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